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Abstract 

Dark matter candidates comprising several sub-states separated by a small mass gap 
Am, and coupled to the Standard Model by (sub-)GeV force carriers, can exhibit non- 
trivial scattering interactions in direct detection experiments. We analyze the secluded 
U(l)5-mediated WIMP scenario, and calculate the elastic and inelastic cross sections for 
multi-component WIMP scattering off nuclei. We find that second-order elastic scattering, 
mediated by virtual excited states, provides strong sensitivity to the parameters of the model 
for a wide range of mass splittings, while for small Am the WIMP excited states have lifetimes 
exceeding the age of the universe, and generically have a fractional relative abundance above 
0.1%. This generates even stronger constraints for Am < 200 keV due to exothermic de- 
excitation events in detectors. 
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1. Introduction 

Many concordant aspects of cosmology and astrophysics present us with compelling evidence 
for a universe in which dark matter (DM) comprises about one quarter of the total energy 
density in the current Universe. However, while we have ample evidence for its gravitational 
interaction, the details of any non-gravitational dynamics if any, and thus its place within 
particle physics, remains obscure. The importance of this question, and the strong motivation 
it implies for physics beyond the Standard Model (SM), has led to an expansive experimental 
program, both terrestrially and in space, aimed at detecting dark matter through non- 
gravitational interactions, namely annihilation, scattering or decay pQ. 

Among the best motivated scenarios for particle dark matter is a generic WIMP - namely 
a weak-scale massive particle, thermally populated during the early universe, and subse- 
quently depleted by a weak-scale annihilation rate [2]. This framework has important impli- 
cations as it suggests the natural scale for self-annihilation and, albeit in a less direct manner, 
gives guidance as to the likely level of scattering with normal matter. Consequently, searches 
for WIMPs present in the galactic halo via their scattering off nuclei in radioactively pure 
underground detectors, have become an integral part of modern subatomic physics pUlUE]- 
Existing searches have primarily focussed on nucleon-WIMP elastic scattering, and the null 
results have placed significant upper limits on the elastic cross-section that are now beginning 
to probe the natural parameter range pQ. 

The chances for successful direct detection depend rather sensitively on whether the 
WIMP is an isolated state, such as real scalar or a Majorana fermion, or a multi-component 
state, such as complex scalar or Dirac fermion. Since the latter allows for the existence of 
vector and/or tensor currents, in many scenarios Z- or 7-mediated scattering can result in 
large elastic cross-sections. However, even a small mass splitting Am between the WIMP 
components - that we will generically denote Xi an d X2 ~ which allows for coupling to these 
currents may significantly alter the scattering signal if Am is in excess of the typical kinetic 
energy of the WIMP-nucleus system. This was first noticed a decade ago [6l [7], where it 
was shown that an otherwise large sneutrino scattering cross-section off nuclei, mediated by 
Z-exchange, could be drastically reduced by splitting the two components of the complex 
scalar sneutrino by Am, with < Am <C m/>. This idea was later exploited in Ref. [8] 
where the inelastic scattering of WIMP states split by Am ~ (9(100) keV was utilized to 
reconcile the annual modulation of energy deposition seen by DAMA with the null results 
of other direct-detection experiments using lighter nuclei. Its of relevance here that the 
signature of inelastic endothermic WIMP scattering, which is enhanced for heavy nuclei [5], 
differs significantly from elastic scattering. 

The possibility of exothermic WIMP-nucleus scattering, i.e. with significant energy re- 
lease, was studied in Ref. [9], where it was shown that a splitting Am ~ O(10 MeV) between 
neutral xi an d electrically charged xf components of the WIMP sector will lead to WIMP- 
nucleus recombination with heavy nuclei via the formation of a bound state between the 
nucleus and the higher-mass negatively charged partner x~2 • A significant amount of energy, 
of 0(1 — 10) MeV, can be released this way via e + , 7, n or v emission depending on the 
particle physics realization of this scenario. Ref. [H] also demonstrated the possibility for an 
O (Am) -1 enhancement in the elastic scattering amplitude due to virtual excitation of the 
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WIMP substructure. 

In general terms, consideration of mult i- component WIMP scenarios with relatively small 
splittings has in the past been motivated on several fronts. On one hand, due to collider 
constraints, supersymmetric scenarios such as the CMSSM generally require an enhanced 
annihilation cross-section in the early universe to avoid over-producing neutralino dark mat- 
ter, and viable parameter ranges usually make use of coannihilation with nearby charged 
states [10]. Another motivation comes from the possibility of the catalysis of nuclear reac- 
tions in the early Universe leading to the resolution of the cosmological lithium problem [11] . 
MeV-scale splittings also allow for the possibility of sourcing the galactic 511 keV line via 
the decay of excited WIMP states [121 [T3J . 

Recently, further motivation for the study of multi-component WIMP scenarios has come 
from claims of a positron excess in cosmic rays above 10 GeV [2], and an enhancement in 
the total electron/positron flux around 800 GeV [IB] . While the interpretation of these 
anomalies remains an active topic of debate, it is tempting to speculate on an origin related 
to dark matter. However, this requires a significant enhancement of the annihilation rate 
associated with cosmological freeze-out, by factors of 10-1000 depending on the WIMP mass. 
This is perhaps most naturally achieved via a new light mediator with GeV or sub-GeV mass 
[T6| [T7] , which allows for a Sommerfeld-type enhancement of annihilation at low velocities in 
the present halo. Decays into these metastable mediators then lead to a dominant leptonic 
branching fraction for kinematic reasons [161 El EE] ■ For the present discussion, the salient 
feature of the enhancement mechanism is that the required Coulomb-like potential arises 
most naturally if the WIMPs are multi-component states, such as complex scalars or Dirac 
fermions, or have at most a small mass splitting between the components [TO1 [T7]F1 

Perhaps the simplest and most natural realization of this proposal [161 E] is via a secluded 
U(l)5 gauge interaction that acts in the dark matter sector and couples to SM particles via 
kinetic mixing with SM hypercharge U(l)y [19]. Such a coupling, kF^ u F^ u , provides one of 
the few renromalizable portals for coupling the SM to new (SM singlet) physics [2D] . The 
particle phenomenology of a (sub-)GeV secluded U(l)5 sector has recently been addressed 
in a number of studies [2T] . with the conclusion that the current sensitivity to k is in the 
range k ~ O(10~ 3 — 10 -2 ). This does not place significant restrictions on WIMP properties, 
as k does not enter the freeze-out cross section [19] , but is interesting in its own right as it 
is close to the natural radiatively generated value [211 122] . 

The question of direct WIMP-nucleus scattering has not yet been addressed in detail for 
U(l)s-mediated scenarios, despite the high level of recent interest. The cross-section was 
calculated to first order in [19] in the limit of small Am, where scattering is mediated by 
a vector current and the WIMP charge radius. This cross-section is generically large, and 
leads to strong limits on k [Ml \TT\, well outside both its natural range due to radiative 
mixing, and consequently the experimental reach [21] . While this 'secluded' regime is of 
interest for indirect detection, it is important to realize as discussed above that such a 

1 We note in passing that pseudo-degenerate xi and X2 WIMP multiplets would also boost the annihi- 
lation cross section by orders of magnitude due to the possibility of (xtx~) resonant scattering, which in 
particular enhances neutralino annihilation into gauge bosons and leptons provided the neutralino-slepton 
system is split by Am ~10 MeV [9]. 
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scattering process can be switched off by even a small mass splitting of the mult i- component 
WIMPs that couple to the U(l)g vector current [THl [32], in complete parallel with earlier 
studies [SI El E] • In this case, elastic WIMP-nucleus scattering may still proceed, but now at 
second order with off-shell excited states. This echoes earlier calculations in [9], where an 
electrically charged partner of the WIMP state is virtually produced and absorbed in the 
elastic scattering process. 

In this paper, we calculate the WIMP-nucleus cross sections for inelastic scattering at 
first order, and elastic scattering at second order, for the minimal U(l)s-mediated model, and 
set constraints on combinations of parameters, such as the mass splitting Am, the kinetic 
mixing parameter k, and the mediator mass my. We also observe that for light mediator 
masses and large mixing, the interaction may become sufficiently strong that perturbation 
theory breaks down. Although Sommerfeld-type enhancement is not important for direct 
scattering with k < 1 and v ~ 10~ 3 in the halo, recombination through the formation of a 
WIMP-nucleus bound state can be kinematically accessible for a relevant range of parameters 
|17j . Another interesting aspect of the minimal U(l)s-mediated WIMP model is the long 
lifetime of \2i the excited WIMP state(s), for small Am < 2m e . For lifetimes in excess of the 
age of the universe, this immediately raises the possibility of collisional de-excitation within 
the detector. Not only are such processes allowed at first order in perturbation theory, but 
they will proceed at a constant rate even in the limit of very small relative velocity. This 
can significantly enhance the energy deposition within the detector, and we use this signal 
to set stringent constraints on the parameter space of the model. These limits employ an 
estimate of the minimal concentration of excited states surviving from the early universe or 
regenerated locally in the galaxy. 

This paper is organized as follows. In the next section we introduce the U(l)s model 
with a mass splitting of the multi-component WIMPs, its low-energy effective Lagrangian 
and the lifetime of the excited states. Section 3 contains the calculation of elastic and 
inelastic scattering cross sections for various parameter choices, and considers the effects of 
WIMP-nucleus binding. Section 4 estimates the fractional cosmological freeze-out abundance 
of excited states, their regeneration in the galaxy, and sets limits on the parameters of the 
model from de-excitation in inelastic exothermic scattering for small Am. We finish with 
some concluding remarks in Section 5. 

2. Multi-component WIMPs with U(l)s-mediators 

The model introduced in Ref. [19] utilizes the kinetic mixing portal to couple the SM and 
(SM singlet) DM sectors via a U(l)g mediator, 

£ = £sm + ^ u (i) s + £dm, (1) 

where £dm describes the WIMP sector and its 11(1)5 interactions. Aj(i) s includes the 
standard U(l) Lagrangian, the kinetic mixing portal and the associated Higgs' sector that 
gives a mass to the U(l)5 gauge boson V^. After this symmetry breaking, the relevant part 
of the low-energy U(l)5 Lagrangian becomes, 

Ai(i) s = -\v^ + + kV^F^, (2) 
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where is the electromagnetic field strength, and a small multiplicative shift by cos 9w is 
absorbed into k. In this paper, for convenience, we choose the Higgs sector to have twice 
the U(l)g charge of the WIMPs. This allows us to introduce two types of (renormalizable) 
mass term in the WIMP action, that itself can be either bosonic or fermionic: 

£q M = ^{iDfjTffj, — m^ip + (XH'tptp + h.c.) fermionic DM, (3) 
4m = {D (D^) - mj0*0 + {\H'<fxf> + h.c.) bosonic DM. (4) 

Here ip (or respectively 0) is the WIMP field in the form of a Dirac fermion (or a charged 
scalar). _D M = d^ + ie'V^ is the usual covariant derivative in terms of the U(l)' gauge coupling 
e', and A is the strength of the dark sector Yukawa interaction. After spontaneous breaking 
of U(l)s, H' develops a vev, and the complex scalar (or respectively Dirac fermion ip), will 
be split into two real (or Majorana) components which we will denote collectively as Xi an d 
Xi with masses rrii = m Xi . The low energy Lagrangian for real scalar WIMPs is then, 

4m = E Qto) 2 - \<xi} + e'V,( Xl d, X 2 ~ X2d,Xi), (5) 

while in the fermionic case, where x% is a t wo component spinor, 

£dm = Yl [^l^AXi - 7,(mXiXi + h.c.) J - ie%(xiVnX2 - X2^Xi)- (6) 

i=l,2 ^ ' 

In both ([2]) and (!5ll6l) . we have suppressed the details of interactions with the physical Higgs' 
as it will not be important for the present discussion. Of more relevance here is that the 
WIMP masses are split by an amount that scales with the Higgs vev v' and the mass of the 
mediator: 

Am = m 2 — mi ~ Xv' ~ —my mi 2, (7) 

e 

where the final inequality reflects the assumed hierarchy, with the mediator much lighter 
than the WIMP(s). For the analysis that follows, this is more relevant than the concrete 
mechanism for producing Am ~ v', as the latter can be achieved within a number of different 
model-dependent realizations. Eqs. (j2J) and fl5|6|) constitute the starting point for the calcu- 
lation of Xi an d X2 scattering off nuclei. However, before proceeding in this direction, we will 
first discuss the requirements on the splitting imposed by the need for Sommerfeld-enhanced 
annihilation in the galaxy, and also the lifetime of the excited WIMP states. 

2.1 Sommerfeld enhanced annihilation 

Sommerfeld enhancement for the scattering and annihilation of WIMPs in the galactic halo, 
with (jra'/v) ^> 1 when »<<a' and small my, originates from the modification of free particle 
wavefunctions due to the impact of the Coulomb interaction when the particle separation 
falls below the de Broglie wavelength (mv)^ 1 . Its then clear that this enhancement will 
only arise (away from possible resonances in the inter-particle potential) if the mass splitting 
between the states is smaller than the Coulomb (or rather Yukawa) potential energy, 

Am < V{r ~ X dB ) =^ Am< ot'mv ~ E kin (—) , (8) 
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where in the latter relations we also assumed that my < mv. Therefore, we find that in order 
to preserve a large enhancement the mass splitting must be smaller than the WIMP kinetic 
energy times the enhancement factor. Thus, for a TeV-scale WIMP in the galactic halo with 
-Skin ~ 1 MeV, requiring an O(10 3 ) enhancement to 'explain' the PAMELA positron excess 
[H] implies the splitting should not exceed 1 GeV. In this paper, we will consider a range of 
mass splittings but generally well below this threshold. 



2.2 Excited state lifetime 



Once produced, for example in the early universe, the longevity of X2 is very important due 
to the possibility for exothermic inelastic scattering in direct detection experiments. The 
decay X2 —> Xi depends sensitively on the size of the splitting. For Am > 2m e , the photon- 
mediated decay to charged particles dominates, and X2 decays on time scales much shorter 
than the age of the Universe for the range in k that we consider. Although this may still 
lead to observable consequences in Big Bang Nucleosynthesis, or during later cosmological 
epochs [23], such effects fall outside the scope of the present paper. The other possibility, 
Am < 2m e , only allows decays to photons and neutrinos within the minimal U(l)g model. 

The X2 ~^ Xi u v decay is mediated by V — Z mixing, and its rate summed over three 
neutrino families is calculated to be, 

4sin 2 ^ GpAm 9 a'n 2 



3157T 3 rriy a 

-53 



3 x 10- GeV x * (*)> ( *5-Y «XY , (9) 



which is much smaller than the inverse age of the Universe, r^ 1 ~ 1.5 x 10~ 42 GeV, unless 
Am > 1 MeV, and/or my < 10 MeV. Note that the rate is suppressed by the ninth power 
of the energy released, as compared to the usual Am 5 scaling for weak decays. This extra 
suppression can be traced back to the fact that V — Z mixing produces a propagator propor- 
tional to q 2 , which is saturated by Am. Similar extra suppression factors are to be expected 
in the decay rate of any SM-singlet WIMPs with derivative interactions to the SM. 

The decay to photons is mediated by the virtual production of an e + e~ pair, leading to 
a loop-suppressed X2 ~^ Xi + 37 decay channel. Unlike fl5j), the decay to photons involves 
V* — > 37 which is phase-space suppressed but does not require a weak transition. Assuming 
a small mass splitting of order 100 keV, we can estimate the rate using the calculated V — * 37 
decay width for on-shell V bosons with my -C 2m e 



^ . , 17 a 3 a'n 2 m% , . 

r v( m„«2m e ) = 1^3W^, (10) 

and replacing my by the momentum q ~ Am of the virtual V. Since the phase space for 
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the 3-photon decay is maximized for momentum q ~ Am, we obtain the estimate, 



r x2-xi+3 7 ~ T v ,^ ?Jl (m* v ~ Am) x a' y—J 

V10- 3 / V 100keV y V / 

which is again much smaller than r^ 1 , unless Am approaches 1 MeV. Although this calcula- 
tion needs to be modified in this limit, we would expect decays above the electron threshold 
to be fairly rapid in any case, as noted above. 

We conclude that for the parameters of interest here, the excited state lifetime can be 
well in excess of the age of the universe. Thus, a primordial excited state population can 
survive from the Big Bang, while scattering in the halo may also re-populate excited states 
in the late universe. Of course, it is important to emphasize that the longevity of the 
excited states in this model follows from the highly suppressed decay rate to neutrinos 
and photons. The presence of additional fields, with a mass below Am, would allow other 
decay channels and thus a significantly enhanced decay rate. For example, one natural 
possibility within the current scenario would be the presence of a very light Higgs' particle, 
so that X2 — ► Xi + h' would be kinematically allowed. However, in this case the longevity of 
Higgs' particle may itself have other cosmological and astrophysical implications that lead 
to independent constraints. 



3. Elastic and inelastic scattering of multi-component WIMPs 

The basic quantity which determines the mass-normalized counting rate in a given detector 
is the differential event rate per unit energy, 

dR AT p x f Vmax , 3 ,, . da 

- Nt ^2L d 3 vvf(v,v E )-—. (12) 



dEf? tyi-y In, . dEft 

Nt is the number of scattering centers per unit detector mass, vn x = vp x /m x is the incident 
WIMP flux, and f(v, ve) is the WIMP velocity distribution in the galactic halo that we take 
to be Maxwellian with escape velocity in the range 500km/s < v E < 600km/s [23]. The 
total number of counts then follows from integrating over the energy bins, and multiplying 
by the effective exposure, e.g. in kg-days, for a given detector [26] . 

The form of the counting rate (112]) isolates the physics of the WIMP-nucleon interaction 
in the differential cross section dajdE R . It has become customary to express the spin- 
independent differential cross section in the form [I], 



da _ m N < ucl 



eff r r rr , r / a ryM 2 



dE R 2v 2 ii 2 n 



f p Z + f n (A-Z) 



nucl 



F 2 (E R ), (13) 



where m N is the nuclear mass, \i n is the WIMP-nucleon reduced mass (which is close to m p 
for electroweak scale WIMPs), and F 2 (Er) denotes a possible recoil-energy dependent form 
factor. This formula presumes that the force mediating the interaction is short-range, so 
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Figure 1: First and second Born amplitudes for Xi-nucleus scattering. 

that the scattering amplitudes on individual nucleons / nuc i can be considered constant. cr^ d 
is then the effective cross section per nucleon, which is now the standard figure of merit for 
comparing the sensitivities of different experiments. 

We now turn specifically to the scattering of quasi-degenerate multi-component WIMPs. 
There are three basic processes that are possible (see Fig. [1]): 

(a) elastic scattering : X(i,2)N — > X(i,2)N. (14) 

(b) endothermic scattering (Q = —Am) : XiN ~~ * X2N. (15) 

(c) exothermic scattering (Q = Am) : X2-N —* XiN- (16) 

The V^-boson vertex in our model always connects Xi an d X2, and thus elastic scattering can 
only occur starting at second order in the Born approximation, while the inelastic processes 
(|T5|) and (fl6|) can occur at leading order, and we will consider these first. 

In order to reduce the parameter space of the model when calculating direct detection 
constraints, we make the well- motivated assumption that the ratio a! /m x is restricted to 
(constant) values which yield the correct relic abundance, i.e. for fermionic WIMPs [T9| [TTj. 

Further details regarding the WIMP relic abundance are discussed in Sec. 4. 



3.1 First-order inelastic scattering 

As discussed in the introduction, the relation between Am and the typical kinetic energy of 
the WIMP-nucleus pair is essential. The first-order Born approximation gives the following 
differential WIMP-nucleus scattering cross section: 

^x 1(2) -x 2(1) _ 4/i ^ aV |k/| F*(q) 



aMf ™ \ki\ (q 2 + m^) 2 ' v ' 

where is the reduced mass of the WIMP-nucleus system, k^/) is the initial (final) mo- 
mentum in the center of mass frame, q is the momentum transfer, and Fjsf(q) is the nuclear 
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charge form factor, — > Z for q^ 1 ^> r^. The magnitude of the momentum transfer de- 
pends on the kinetic energy in the c.o.m. frame E kin = kf/(2/i7v) = k 2 /(2fi N ), the scattering 
angle 9 and the mass splitting Am, 



^ = 2E kin (l-Jl T ^ C os6 ) =F Am, (19) 

where the minus sign corresponds to Xi ~ > X2 and the plus sign to X2 Xi scattering. 
The recoil energy of the nucleus in ([T2]) is then Er = q L 2 /2m^. The relation ffT9]) illustrates 
that for the endothermic process ( [15]) . = Am, while for the exothermic scattering ( fT6l) . 
£™ n = where q 2 (E kin = 0) = 2Am/j N . I n the limit Am < £ kin , Eq. ([19]) reduces to the 
standard elastic scattering relation, q = a/2P(1 — cos#). In the limit E^fi^ <C my and 
neglecting the form factor dependence, ^ {E^hn) 1 ^ \ the total cross section reduces to 
the elementary formula [19] : 

16nZ 2 aa'K 2 fi 2 N 

°"xi(2)^x2(i) - m ^ ' y Z(} ) 

which can be interpreted as the scattering of a WIMP due to its finite electromagnetic 
charge radius [27]. In the opposite limit of small m v , E kin fi N 3> m v , the scattering becomes 
Rutherford-like, and the total cross-section saturates at the minimal momentum transfer, a ~ 
which in an experimental setting would be determined by the lowest detectable recoil 
energy Er. The generalization of ([20]) to a finite mass splitting in the limit E kin fj,^, Am/iN "C 
my is straightforward: 



16nZ 2 aa , K 2 fi 2 N / Am 

a w*«» = < V T e£ { ] 

For slow incoming particles with Am ^> E^, we see that for exothermic scattering the 
factor y/^~- oc 1/v so that the rate is enhanced by the familiar a ~ 1/v dependence of the 
cross section. 

To set constraints on the parameters of the U(l)s model, we insert the differential cross 
section ffTB"]) into Eq. ([12]) and calculate the 90% Poisson confidence limits (CL) from direct- 
detection data. In order to facilitate as simple a comparison as possible between constraints 
from endothermic, exothermic, and elastic scattering, in this paper we will restrict ourselves 
to data from the CDMS experiment [3] . Conservative limits result from the assumption that 
the galactic dark matter component is dominated by \ii an d so we focus for the moment on 
the endothermic first-order process. We shall return to constraints following from exothermic 
X2 Xi scattering in the next section. As can be seen in Fig. [2J the strong sensitivity to 
K,/my at Am = noted in [T6] [T7] diminishes as Am increases. Moreover, for a given recoil 
energy, there exists a certain critical value of mass splitting above which kinematics forbids 
the endothermic inelastic scattering. Thus if Am is large enough the sensitivity is completely 
lost, and this occurs for Am > 190 keV for the parameters used in Fig. [21 

3.2 Second-order elastic scattering 

With a larger mass splitting Am > MeV between the WIMP xi an d its excited partner 
X2, kinematics forbids the first-order scattering process xiN — * X2N . However, elastic 
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m x =200 GeV 
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Figure 2: Endothermic inelastic scattering constraints. 90% CDMS confidence limits on k 
as a function of vector mass my for a 200 GeV WIMP, with a' chosen to yield the correct 
thermal relic abundance. We show constraints from the endothermic scattering XiN ~^ X2N 
for Am = (solid), 100 keV (dashed), and 150 keV (dotted). No endothermic events are 
expected in this case for splittings above Am ~ 190 keV. 

scattering (|T4l) can still proceed via double V — 7 exchange as depicted in Fig. [TJ For a 
light U(l)s vector mediator with my < 50 MeV or so, the WIMP-nucleus interaction range 
is comparable to or larger than the typical nuclear radius tn, leading to saturation of the 
virtual momenta at | q virt | ~ my < r^ 1 and resulting in an (9(Z 2 )-scaling of the WIMP- 
nucleus scattering amplitude. This implies that such an amplitude cannot be decomposed 
into a sum of scattering processes on individual nucleons. Consequently, in this regime 
Eq. (|T3j) is applicable only as a definition of the effective nucleon amplitudes that would lead 
to the equivalent WIMP-nucleus cross section. However, once the range of the mediating 
force falls below 1 fm and my starts approaching a GeV, WIMP scattering starts probing 
the internal nuclear structure, and ultimately the scattering amplitude does become a sum 
of amplitudes for scattering on individual nucleons, as in (1T31) . Further increase in my above 
the GeV scale allows the U(l)s sector to be integrated out leading to a WIMP-quark effective 
Lagrangian, that can be used to calculate f p ( n )- 

In what follows, we calculate the second order elastic scattering amplitude, keeping m v 
under a GeV. For light and intermediate- mass U(l)s mediators the WIMP-nucleus cross sec- 
tion can be straightforwardly calculated using a non-relativistic potential scattering treat- 
ment in the presence of a nuclear charge form factor. For m v x well below 1 fm, analo- 
gous formulae can easily be derived for WIMP-nucleon scattering. The wavefunction of the 
WIMP-nucleus system under the influence of a central potential is ip(x.) = e* ki ' x +/(#)e* k/ ' x / r, 
where f(8) is the scattering amplitude. The second-order Born amplitude for scattering from 
a potential ^(x) can be written in general as 

f(Q) = -0± /d 3 xd 3 x'e^ x K(x)G'(x-x / )V(x')e iklX ', (22) 

27T J 

where G = (E^ — E' v )~ l is the nonrelativistic propagator for the ^-nucleus system, which 
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in momentum space becomes, 



G(p) = (ym^ + k 2 + yjm% + k 2 ) - ^mj s + p 2 + y/m% + p 2 J 



k 2 p 2 

—Am + 



2//jv 



(23) 



where p is the intermediate momentum. For a large enough splitting we can neglect the 
kinetic energy terms in the propagator, and this occurs for Am > E^ n ~ 100 keV for a 
typical nucleus with mass of order = 50 GeV. Retaining only the constant (Am) -1 term 
in the propagator allows us to use the completeness relation and remove the integration over 
intermediate states. Thus the amplitude greatly simplifies and is expressed via the square 
of central potential V(x) = V(r): 

2 /^V ? „ T r2, 



f(q) = —i — / drrV (r)smqr. (24) 
qAm 



where q = a/2A; 2 (1 — cos 9) and 9 are the momentum transfer and the scattering angle in the 
c.o.m. frame, as before. The expression (1211) is equivalent to the first-order Born formula for 
the effective potential V e g = V 2 (r)/Am. 

The potential can be written in general as 



d 3 p F(p) 



V(r) = KVaa' \ ~ ~ ■ 2 e 
(Ztc) 6 p z + m 



ipx 

V 



KZ^fart f 00 (F(p)/Z)p sin pr 
= ~^h~J dP v 2 + m\ {2b) 

where F(p)/Z is the charge-normalized nuclear form factor, with F(0)/Z = 1. To compute 
this potential, we must specify F(p). Using a model with a uniform charge distribution 
within a sphere of finite radius radius r N , gives 

F(p) = mWN \ (26) 
pr N 

where ji is a spherical Bessel function. This simple model of the nucleus allows us to 
analytically compute the potential in Eq. fl^Bll 2 !: 

3ZkV^ (/(^^)l,.< r . JV = l-(l + ^)e" :c ^^. 

v ( r ) = — 3 — f(m v r N ,m v r) where < r e _ f (27) 

m V r N [f( X ^)\r>r N = ( X COshs — SUlll x)— —. 



2 Note that for the first-order endothermic and exothermic inelastic processes we utilize the more common 
Helm form factor 28J, which contains an additional exponential factor exp(— s 2 p 2 ) relative to Eq. (l2Gp . 
The parameter s is the nuclear skin depth with a typical value of 1 fm, corresponding to energy scales of 
order 200 MeV. For the case of coherent elastic WIMP-nucleus scattering considered here, the interaction 
range is approximately rriy 1 ~ > s, so that the simple solid sphere model of the nucleus is a reasonable 
approximation. 
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This is essentially the 'Yukawa potential' of a hard sphere, and it is clear from (1271) that 
in the limit m v — > at finite r N , the potential reduces to the usual electrostatic potential 
created by a spherical charge distribution, while for r N — > 0, V(r) ~ exp(— m v r)/r. Taking 
both limits simultaneously reduces fl27j) to an ordinary Coulomb potential, in which case the 
scattering amplitude (jUj) simplifies dramatically to give, 



_ „ u = ^ff^)\ j^_\ t (28) 

q/\m k z {/\m) z \r^my J 

where the final integral over momentum transfer is approximated by a logarithm given the 
assumption my <C q <C r^ 1 . Because the range of validity for (1281) is almost non-existent, 
and the nuclear form factor has been neglected, we will not make use of it for deriving 
constraints on the model. However, it has the merit of illustrating the Z 4 -scaling of the 
cross section that translates into the following Z-independent prediction for the effective 
cross section per nucleon 

-49 fa'V ( k \4 (\ MeV\ 2 „ 

a nucl ^10- 42 cmx f-1 {-^^j for small my (29) 

where we have inserted typical values for various parameters such as Z/A ~ 0.5, the char- 
acteristic recoil energy, etc. We can see that direct detection experiments can potentially 
have strong sensitivity to k for light mediators, due to the resulting long range interac- 
tion. Eq. (1291) represents the sensitivity in the 'best case', which erodes rather rapidly with 
increasing my. 

To derive actual limits on the parameters of the model, we write the differential cross 
section as 

±- = ?=^|/( Bj 0|». (30) 
dE R /i> 2 

Inserting Eq. (|30|) with full form factor dependence into (|T2|) . we calculate the constraints 
on the parameter space using the results from the CDMS experiment. These constraints 
are shown in Fig. [3] for two cases of my = 10 and 100 MeV. It is easy to see the reduced 
sensitivity to k as My 1 becomes smaller than the nuclear radius. 

For larger mediator masses, one can easily observe this change in the Z-scaling. Once 
rriyr N ^> 1, then inside the nucleus the potential (TSTJ becomes V(r) ~ ZZk\J aa'm v 2 rjf 3 , 
and the amplitude at small q has the following scaling with Z: f(q = 0) ~ Z 2 r^ 3 ~ Z. 
The switch from Z 2 to Z in the amplitude signifies the loss of coherent nuclear response, 
and at that point the process is best described by scattering on individual nucleons. Taking 
my ~ GeV, we calculate the WIMP-nucleon elastic scattering amplitude, following the same 
procedure as developed for the nucleus. We use a proton form factor F p {p) = (1 + a 2 p 2 )~ 2 
with a ~ 0.84 fm related to the charge radius of the proton. The amplitude in Eq. ( 1241) . 
the cross section and the rate can all be straightforwardly computed, and we obtain the 
estimate: 

42 {u'Y ( k \ 4 flO MeV\ 2 „ 

a nucl ^10- 42 cmx (-J (— ) [~^ r ) for my ~1 GeV. (31) 
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Figure 3: Elastic scattering constraints: 90% CDMS confidence limits on k as a function 
of WIMP mass m x for a mass splitting Am = 10 MeV, with a' chosen to yield the correct 
thermal relic abundance. We show constraints from the elastic scattering XiN — > XiN for 
m v = 100 MeV (solid) and 10 MeV (dashed). 

It is apparent that the loss of coherence in the elastic scattering process means that the 
sensitivity to k in this case is no better than particle physics probes with SM particles. The 
dependence on my of the sensitivity is exhibited in Fig. @J showing the transition around 
m v ~ 1/rjv- 

Finally, for a fermionic WIMP, a loop-induced effective coupling of the nucleon spin to 
an off-shell photon become possible. This can be interpreted as an effective anapole moment 
of the WIMP [27J IS]- Although such a coupling will be proportional to the first power of k, 
it has no implications for the direct detection, primarily because the scattering amplitude is 
spin-dependent and this not coherently enhanced, as well as being proportional to the small 
relative velocity of the WIMP-nucleus system. 

3.3 Non-perturbative scattering and WIMP-nucleus binding 

In the previous section we assumed that the Born approximation is applicable, which may 
not always be the case if my is small and k is taken to be large. For small Am and my, the 
criterion for the applicability of the Born approximation is Zy/aa'n/v <C 1. Since Za could 
maximally be of order one, we see that for galactic WIMPs with v ~ 10~ 3 , this is satisfied 
for any k < 10~ 3 . We note that for small my the exclusion contours in Fig. H] are well 
below this threshold and thus well within the range of validity of the Born approximation. 
Nonetheless, it is intriguing to consider whether V — 7 exchange may lead to the formation 
of a bound state as suggested in Ref. [TTj. 

For a suitable choice of parameters, it is indeed possible for the Yukawa-like potential 
between the WIMP and a large nucleus to support a bound state, or have a quasi- stationary 
state just above the continuum threshold. The potential arises from V — 7 exchange, and in 
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Figure 4: Elastic scattering constraints: 90% CDMS confidence limits on k as a function of 
vector mass my for a 200 GeV WIMP with mass splitting Am = 10 MeV, where a' is chosen 
to yield the correct thermal relic abundance. We show constraints from elastic scattering 
X\N —>■ xiN , where the solid line interpolates between the nuclear and nucleon scattering 
descriptions. 



the limit my <C 1/r/v takes the form 



.Vac? 



V(r) = -KZ^—e- m "\ (32) 

In the limit of small Am, the existence of at least one bound state then requires (for <C 

m x ) 

K > 1.4— , (33) 

which can reach interesting parts of the parameter space for light mediators, and large nuclei, 
e.g. for the Germanium target at CDMS, a bound state is possible for k > 10 -3 for a' = a 
and my ~ 10 MeV. 

The binding energy is rather small, varying from the Coulomb limit Eb = (K,Z) 2 aa'm x /4:, 
which could be up to 0(100 keV) for large (kZ) ~ 0.1 down to Eb = at the threshold fl33|) . 
For the situation where the Dirac WIMP is split by Am into Majorana components, this 
potential will inhabit the off-diagonal terms in a 2 x 2 potential matrix, thus for a mass 
splitting Am ~ 100 keV, its apparent that only large nuclei will allow for recombination. 

If kinematically accessible, since the recombination rate for \N — > (xN) + 7 is large, 
we can conservatively regard the threshold (1331) as a constraint on the parameters of the 
model for small splittings. This would most likely arise through the formation of too large 
a relative abundance of anomalous heavy elements, but it would also produce a distinctive 
signature in underground detectors. Nonetheless, it is worth remarking that that the strong 
dependence on the atomic mass means that the most stringent constraints on the abundance 
of anomalous heavy isotopes of light elements will, as was the case for the pseudo-degenerate 
scenarios of [9], be less important. Finally, if Am is increased or k decreased, the bound state 
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can be pushed into the continuum, and may lead to a resonant feature in WIMP-nucleus 
scattering. However, this again will likely be well inside the region excluded by Fig. HI 



4. Inelastic scattering and de-excitation constraints 

The presence of an ambient population of excited \2 states in the galactic halo would allow 
for exothermic inelastic \2 Xi scattering (fl6l) in detectors. Given the enhanced rate, this 
can potentially act as a novel and powerful probe of multi-component WIMP scenarios in 
which a sizable excited state population is present. However, these constraints will depend 
sensitively on Am, as a large O(MeV) splitting would result in too great an energy deposition 
within the detector, E ~ Amm x / (m x +m7v), falling outside the fiducial recoil energy window 
for most searches. Moreover, larger values for Am may considerably shorten the lifetime of 
X2, as in Eq. (jTTI) . reducing the local number density. Therefore, in this section, we will 
concentrate on a relatively small splitting: Am < 500 keV. 

While the overall abundance of \2 arid xi is fixed as soon as chemical equilibrium is 
lost in the early universe, namely after the temperature drops below 0.05m x , the relative 
abundance of X2 arid xi can vary. The possibility of exothermic scattering in detectors 
depends on the ambient number density n 2 of the excited X2 population in the galactic halo, 
which arises from three main sources, cosmological, galactic, and local: 



The cosmological abundance is regulated by the X2 lifetime and freeze-out of the X2 — » Xi 
rate in the early universe, while the galactic source is related to the inverse possibility for 
Xi ~ ¥ X2 up-scattering in the galaxy, with energy supplied by the WIMP kinetic energy 
[T2| [T3] . The local source could originate from the scattering of xi on heavy nuclei in the 
Earth's interior, e.g. XiPb — > X2Pb. It is relatively insensitive to the X2 lifetime, provided 
that it is longer than the WIMP time-of-flight through the Earth. However, it is easy to 
see that the weak scale WIMP-nucleus cross section cannot lead to more than an O(10~ 8 ) 
excitation probability. Thus the cosmological and galactic sources are in general far more 
important, but depend crucially on the size of the X2 Xi inter-conversion rate, which may 
in general have a number of contributions: 

• Double (de-) excitation: 



This is the most generic possibility, and the only one which will be important in the 
minimal U(l)s model for a generic range of parameters. Note that it depends sensitively 
on the number density of X2 states, and thus will tend to freeze out - at the latest - 
when the temperature drops below Am. This will allow us to determine the minimal 
fractional abundance of X2 relative to Xi in the 11(1)5 scenario. 

• SM thermalization: 



(c) , (<?) , (0 
n,2 = n 2 + n 2 + n-2 ■ 



(34) 



X2 + X2 <-> Xi + Xi- 



(35) 



X2 + SM <-> xi + SM. 



(36) 
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For this process to compete with (I35I) requires a significant interaction rate between the 
WIMP and SM sectors^, and amounts to maintaining thermal contact of the WIMPs 
with the SM bath down to low temperatures T 7 < Am. Since we take Am < m e , this 
may only happen if the rate for scattering off electrons is larger than the weak rate, 
An n\/ aa' / m v > Gf, because for T 7 < m e the number density of electrons decreases 
exponentially However, for the fiducial range of values that we consider, k\J aa'/m v 
is broadly comparable to Gj?, and the scattering off electrons will have the effect of 
maintaining thermal equilibrium between the WIMP and SM sectors [U [2H] , possibly 
all the way to T 7 ~ MeV, but will not in itself lead to depletion of the \2 abundance. 

• Dark thermalization: 

X2 + X^ X i + X. (37) 

This process assumes a more complex dark sector with additional light degrees of 
freedom. Given that the dark sector decouples relatively early from the SM bath, 
and subsequently cools more quickly, if this additional de-excitation rate can remain 
in equilibrium well below this decoupling scale then it could allow for a significant 
depletion of the excited state fraction. Such a process then becomes intertwined with 
the possibility for X2 — ► Xi + X decay. This is not relevant for the minimal U(l)s 
scenario for my > 0(MeV) as considered here, and so we will not consider it in 
detail. However, it is worth noting that the presence of additional light states may 
have further implications for BBN and cosmology, particularly if this interaction rate 
is large as would be needed for this process to be important. 

In the minimal U(l)s scenario, for Am < m e , the double de-excitation process fl35l) is 
generally the most relevant and we will now consider the freeze-out in more detail, in order 
to estimate the minimal fractional abundance of X2 relevant for direct detection. 

4.1 Estimate of the fractional abundance of excited WIMPs 

In the early universe, chemical freeze-out occurs at relatively high scales, where the distinc- 
tion between xi an d X2 is unimportant. Since we have m x 3> my, the primary annihilation 
process involves XX ~^ W where the two V bosons are on-shell and subsequently decay to 
the SM via kinetic mixing with the photon. For my <C m x we have [TSJ [T7| 

Mann = ~> 2.4 X KT^CHlV 1 , (38) 

where the latter relation follows from ensuring that the WIMPs have a relic density that 
saturates the measured value of VL DM . For our purposes, chemical freeze-out for xi an d X2 
will occur at a high temperature scale Tf ~ m x /20, and so we can simply take (138]) as a 
constraint relating the coupling a' to m x as in (fT7|) . 

If the X2 lifetime exceeds ru, as discussed in Sect. 2.2, and the rate for inter-conversion 
( 1351) is slow - dropping below the Hubble expansion rate before the average WIMP energy 

3 We are grateful to D. Morrissey for pointing out the relevance of SM thermalization. 
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falls below Am - the inevitable prediction is n 2 /n\ ~ 1. For example, if in the U(l)s model 
the range of the mediator m v x is very short (e.g. weak scale), X2 and xi wm be equally 
abundant. However, given a relatively long interaction range for small my, the reaction (1351) 
could remain in equilibrium down to WIMP energies of E < Am, resulting in a significant 
depletion of X2 states. 

After chemical decoupling, the WIMPs remain in thermal equilibrium down to lower 
temperatures, scattering off U(l)5 vectors in the dark sector via the Thomson-like process 
XV ~^ XV and/or off SM charged particles. Scattering off V^'s becomes inefficient once the 
temperature drops below my, and most of the vector particles decay. The ^-dependent 
scattering of WIMPs on SM charged particles can be straightforwardly computed e.g. in 
the temperature interval 1 MeV < T < 100 MeV where scattering off electrons is dominant, 
with a rate scaling as K 2 aa'rriy 4 T 5 . After these processes fall below the Hubble rate, the 
WIMPs are no longer in kinetic equilibrium with the SM thermal bath and begin to cool 
more rapidly, maintaining a quasi-thermal spectrum with temperature T x : 

T X = T 7 2 /T, for T 7 <T*. (39) 

For very small k, scattering on V dominates, and T* ~ my. However, if the rate of re- 
scattering on electrons is comparable to the weak rate, i.e. A-KK^aa'Y^my 2 ~ 10 -5 m,p 2 , 
decoupling is postponed to a temperature T* ~ MeV. 

Once kinetic decoupling from the SM bath is complete, the inter-conversion process ( l35l) 
is capable of driving down the fractional abundance of X2, provided that it is faster than the 
Hubble rate, and the energy of the emerging non-thermal xi s is quickly re-distributed in the 
WIMP sector. If such a quasi-thermal state is maintained below T x < Am, the fractional 
abundance of excited states will be exponentially suppressed, 

n 2 ( Am\ , . 

~ exp -— - . (40) 



F I rri I ' 

Due to the finite rate for the X2X2 X1X1 de-excitation process, the approximate freeze-out 
condition is, 

H(Tf) = [n 2 (a 22 ^ n v)} T f, (41) 

where the SM and dark sector freeze-out temperatures are related via (1391) . This results in 
the following estimate for the fractional freeze-out abundance: 

n 2 _ ^-1 Jh_ H ( T J) (42) 



n 2 + ni Vt DM m p [n 1 (a 2 2^nv)\ T f 

where rji, = 6.2 x 10~ 10 is the baryon-to-photon ratio, Qb/^DM — 0.2 is the ratio of baryonic 
to DM energy densities, m p is the proton mass, and n^(T) = 0.24T 3 is the photon number 
density. The Hubble rate is given by HiT) ~ l.lg^M^T 2 , while M P1 = 1.2 x 10 19 GeV and 

g* ~ 10. 

The minimal value for n 2 jn\ is achieved when the Xi X2 inter-conversion cross section 
is maximized. Since we are working in the regime a'/v > 1 with small my, perturbation 
theory is not applicable, and the Schrodinger equation should be solved numerically as in 
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the recent paper [30] . However, for our discussion it suffices to saturate cr 22 ^ii by the s-wave 
unitarity limit or by the range of the force carrier, m y 2 , whichever is smaller: 



<«n ~ ^ for k>m v 



{a 2 2-*nv) < 



a/a 



3/2rr 

m x ' T 7 



(43) 



Combining (j4"2"l) and (|43p . we obtain the minimal fractional \2 abundance surviving from the 
early cosmological epoch: 



n2 



(c) 



10 



-2 



300 GeV7 



5/2 



10 MeV\ 1/2 



T., 



(44) 



where the abundance is also implicitly bounded from above by [^2,Mi]max = 1- Note that 
if thermalization via electron scattering is negligible, then in order to maximize the rate 
one can chose my ~ k ~ (in^Am) 1 / 2 at WIMP energies comparable to Am, so that the 
optimal value for T* in (1441) is T* ~ my ~ (m x Am) 1 ^ 2 . More generally, the estimate (jUJ) 
indicates that for a TeV-scale WIMP, the fractional abundance of X2 does not drop below 
5%, independent of the strength of the de-excitation reaction. We should note that more 
accurate computations of the excited fraction n 2 /ni are certainly possible, but the estimate 
will suffice to obtain an estimate of the direct detection rate. 



In addition to this cosmological source, there is also the probability of endothermic up- 
scattering in the galaxy, adding to the fractional abundance in ( 14"4"|) . Since the rates for the 
forward and backward reactions are related, we can estimate the excited state abundance 
from galactic scattering to be, 



Tl2 



Cff) 



r 'mt 



m v 



(o-11^22w)gal, 



(45) 



where p x ~ 0.3 GeV cm -3 is the local dark matter energy density, v ~ 10~ 3 , and r int the 
'integration time', equal at least to the age of the Milky Way, 13 bn yr, or the lifetime of X2, 
whichever is smaller. Assuming a lifetime in excess of 10 bn yr, we obtain a simple estimate, 



U2 

ni 



i(g) 



~ 10~ 4 x 



300 GeV 



x exp(-Am/T eff ), 



(46) 



where T e g is the effective WIMP temperature in the galactic halo. This suggests that the 
galactic source is somewhat subdominant, but can in principle compete with the cosmological 
source if the WIMP mass is ~ 100 GeV. 



4.2 Constraints from direct detection 

Making use of Eq. (jUJ) to determine the minimal fractional \2 abundance, and the results of 
the previous section, we are now able to calculate the exothermic scattering rate as a function 
of the parameters in the U(l)s model. To be specific, we choose my = 1 GeV, and vary Am 
and k, recalling that elastic scattering processes do not place significant constraints on k for 
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Figure 5: Exothermic inelastic scattering constraints. 90% CDMS confidence limits on k 
as a function of mass splitting Am for a vector mass my = 1 GeV, where a' is chosen to 
yield the correct thermal relic abundance. We show constraints from exothermic inelastic 
scattering X2N Xi-N for WIMPs with masses m x = 100 GeV (solid) and 1 TeV (dashed). 
The constraints rapidly deteriorate for large Am as in this case most scatterings will have a 
large nuclear recoil well above the maximum detector sensitivity of Er = 100 keV. 

such large values of my. The CDMS constraints for a 100 GeV and a 1 TeV WIMP are 
presented in Fig. [51 It is clear that for small values of Am exothermic scattering places very 
stringent bounds on values of the mixing parameter k. Also, as one expects, the constraints 
start to deteriorate for splittings Am > 100 keV, as these exothermic events will generally 
have a large recoil energy, En > 100 keV, outside the window of most direct detection 
experiments. 

Another intriguing feature of Fig. is the increased sensitivity of direct-detection experi- 
ments for heavy WIMPs. Naively the differential rate (Tl2|) falls as m~ , but in our case there 
is a hidden m x dependence in the sensitivity coming from both the necessity of larger U(l)s 
couplings a' to obtain the appropriate cosmological DM abundance, and the WIMP mass 
dependence contained in our estimate of the minimal fractional abundance in Eq. f j44l) . 

As an application, we consider the feasibility of the minimal U(l)s scenario to account for 
the DAMA anomaly [5] in light of constraints coming from both endothermic and exothermic 
scattering. The inelastic dark matter scenario [8] was proposed to reconcile the null results 
of direct detection experiments with the annual modulation signal observed by the DAMA 
experiment, and multi-component WIMPs in the U(l)s scenario provide one realization of 
this idea. To calculate the DAMA preferred regions, we utilize a \ 2 goodness-of-fit test [31 J, 
including data from the first 12 bins between 2 keVee and 8 keVee [5] . We show the combined 
constraints from endo- and exo-thermic processes from CDMS in Fig. O for two values of 
the WIMP mass, 100 GeV and 1 TeV, along with the DAMA preferred regions at the 90% 
and 99% CL. Since the relevant range of n, my and a' considered in this paper implies 
a rate for WIMP scattering off the thermal electron-positron bath which is close to weak- 
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Figure 6: Constraints from endothermic (solid) and exothermic (dashed) scattering in the 
CDMS detector at 90% confidence level for a 100 GeV WIMP and 1 GeV U(l) s vector, with 
a' chosen to yield the correct thermal relic abundance. We show constraints for a 100 GeV 
WIMP (left panel) and 1 TeV WIMP (right panel). The complementarity of both constraints 
is clearly seen. The light (dark) shaded region corresponds to the DAMA 99% (90%) CL 
preferred region, which is excluded by both constraints. 

scale even for k ~ 10~ 4 as discussed above, the relevant choice of decoupling temperature is 
T, ~ 0(1 - lOMeV). 

One observes that in various parameter ranges exothermic scattering can in principle 
provide more stringent bounds on the parameters of the model than endothermic scattering 
or second-order elastic scattering for small Am, as shown in Figs. [2]and|H One should keep 
in mind, however, the model-dependence of such a conclusion as discussed earlier in this 
section. Indeed, any extension of the minimal U(l)s model with additional light particles that 
could facilitate the decay/de-excitation of \2 may result in the absence of any significant \2 
population. The elastic and endothermic scattering constraints are, in contrast, completely 
independent of details concerning the \2 lifetime and/or de-excitation rate. 

Finally, it is also worth noting that another possible signature is via ionization due to 
exothermic WIMP scattering on atomic electrons in the detector, which will be essentially 
monochromatic as it is practically independent of the WIMP velocity, with a cross section 
that scales linearly with Z. 

5. Discussion and conclusions 

In this paper, we have carried out a systematic study of the nuclear scattering of multi- 
component WIMPs forming part of a SM singlet dark sector, focussing on its minimal im- 
plementation via a U(l)s mediator. The mass splitting of the components by Am leads to 
rich structure for scattering phenomenology. While forbidding first-order elastic scattering, 
the second-order Born cross-section still implies significant sensitivity for direct detection 
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experiments such as CDMS and XENON for a wide range of mass splittings. Moreover, 
for small splittings of (9(100) keV, first-order inelastic scattering provides far more stringent 
constraints, particularly as in this regime the excited states may have a lifetime exceeding 
the age of the universe and thus a residual population in the halo can allow exothermic 
down-scattering in the detector. This exothermic process leads to additional constraints on 
specific 'inelastic DM' scenarios, via scattering on relatively light target nuclei that would 
otherwise have no sensitivity for larger Am. 

In rather general terms, it seems there is a clear tension between the constraints discussed 
here, associated with enhanced scattering rates, and the presence of relic populations of 
excited states that may be relevant for other aspects of dark matter phenomenology. We 
will conclude by discussing some specific implications in this vein, but going beyond the 
specific U(l)s scenario: 

• Inelastic DM and DAMA: It was noticed a few years ago [8], that a mult i- component 
WIMP scenario with a small 100 keV splitting could resolve the tension between the 
annual modulation signal observed by DAMA [5] , with the fact that the required cross 
sections were seemingly ruled out by the null results of other experiments with lighter 
materials. Despite further null results from XENON [1], recent analyses suggest there 
is seemingly a small parameter range where all experiments can be reconciled [32J . Our 
simple observation here is that with such a small splitting, the lifetime of a SM singlet 
WIMP excited state will generically be long, and thus the possibility of exothermic 
scattering of a residual excited state population will be an important constraint as in 
Fig. [H]for the U(l)s scenario. Ameliorating this constraint would require consideration 
of models with further decay channels, e.g. via adding a SM charge to the WIMP (as 
for a sneutrino), or via additional light states. 

• Multi- component states and the INTEGRAL 511 keV line: Attempts to attribute the 
well- measured 511 keV line from the galactic center to WIMP interactions generically 
involve metastable states with decays to positrons. This may be via delayed decays of 
a relic |13j . or via excitation and decay as in the 'exciting dark matter' proposal |12j . 
In the latter case, utilizing the WIMP kinetic energy in the halo requires a scattering 
cross section well in excess of s-wave unitarity [321 ED] unless optimistic assumptions are 
made concerning the galactic halo profile. Recently, it was suggested that introducing 
a three-component WIMP sector {xi?X2,X3} [3H] would resolve this issue by utilizing 
a large relic population of the middle state X2, with a small 100 keV splitting from the 
highest state X3- Scattering could then excite this transition with a reasonable cross- 
section, with the subsequent rapid decay to the ground state Xi liberating 0(1 MeV) 
and sourcing the 511 keV line. In this case its clear that the lifetime and fractional 
abundance for xi must be large, suggesting an analogous problem with the large cross 
section for exothermic scattering in direct detection. However, it is worth noting that 
in this scenario the energy release would be of (9(1 MeV), which may in fact be too 
large to be easily identifiable given the current fiducial recoil energy range. 

• Direct detection of exothermic scattering: These issues raise the important question of 
the detectability of exothermic scattering in direct detection experiments. It is clear 
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that for larger splittings, Am > 300 keV, the peak of the recoil energy spectrum tends 
to fall outside the conventional fiducial range. It would clearly be of interest to know 
whether existing experiments are able to expand the conventional 10-100 keV fiducial 
region, in order to provide more stringent constraints on these scenarios. 

Note Added: While this work was being completed, the preprint [33] appeared on the 
arXiv, which also considers exothermic scattering through de-excitation, and thus has some 
overlap with the discussion in section 4. 
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